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1. Introduction 


Albania is a small country of only 28 787 km? surface area and 
around 4500000 inhabitants, situated in the southwest part of 
the Balkan Peninsula. Like the other Balkan countries, Albania 
is located next to the subduction boundary between the African 
plate and the Euro-Asiatic plate. This setting makes the pres- 
ence of geothermal resources possible. Surface manifestations 
of geothermal resources are found throughout Albania, ranging 
from the region of Peshkopia in the northeast, where hot springs 
with water temperature of about 43°C and inflow above 141/s 
are found, through the central part of the country with different 
sources (including the springs of Llixha-Elbasan) with temper- 
atures above 66°C, to the Peri-Adriatic depression (see Fig. 1) 
with a number of wells (drilled for oil and gas research) pro- 
ducing water with temperatures around 40°C at variable yields. 
The thermal waters in Albania are only used for balneology. 
This form of use dates back to early in history, or to the time 
of the Roman Empire (i.e. the Sarandaporo’s thermal baths). So 
far, the geothermal resources have not been utilized for other 
purposes, such as space heating. Estimated temperature mea- 
surements based on different geothermometers indicates that the 
temperature of the waters in the formation of the Llixha reser- 
voir may be above 220°C. The reservoir is believed to be in the 
depth interval of 4500-5000 m. Tirana’s geothermal heat sources 
are: 


e Underground waters of the Tirana quaternary depression; 
e Underground waters of the tortonian’s sandstones; 

e Heat of the peri-superficial quaternary or tortonian layers; 
e Limestone’s and dolomites saturated with artesian waters. 


The following is addressed in the material: 


e The general geological conditions in the areas; 

e A review of the theoretical basis of heat transfer; 

e Finite-volume modeling of the whole geothermal system down 
to 5000 m depth, incorporating both thermal convection and con- 
duction, based on a simple boundary conceptual model; 

e House heating calculations; 

è Borehole heat exchangers calculations; 

e Economical analyses. 


All this aims at demonstrating that the thermal water flow- 
ing from the Llixha springs is usable for direct utilization and the 
borehole heat exchanger heating systems are economically feasi- 
ble, despite their elevated preliminary cost. This utilization would 
mitigate the electricity supply to the region, help improve living 
conditions for the local community and protect the environment 
of the region. 


2. Geological background of the Llixha-Elbasan and Tirana 
region 


2.1. The geological structure of Llixha region 


The Llixha region is situated southwest of Elbasan. The region 
is well known for its thermal springs, appreciated since ancient 
times for their curative properties. A geological map of the 
region and a representative cross-section (I-I) are presented in 
Figs. 2 and 3. 

The region under study is south of Shkumbini river valley. 
The surface relief increases rapidly up to intermediate eleva- 
tion (300-500 m). In the western part of the region, a system of 
hills declines gradually in the Cérriku field. The region has a rich 


hydro system of small streams and many underground water- 
systems. The flow rate for the underground waters varies from 
100 to 2001/h in Thanë up to 50001/h in Tregan. Generally, the 
formations are composed of flysch with a diverse and chaotic mor- 
phology. This is an inhabited area, with small villages clustered 
around the thermal waters no more than 2-3 km from each other 
[1]. The region lies between two tectonic regions; the transver- 
sal Vloré-Elbasan-Dibér and the longitudinal Leskovik-Drini river 
bay. Both of these connect the lower part in connection with 
the northern and western part of the country. In the context of 
Albanian tectonics, the region represents the western part of the 
Kruja tectonic zone. The formation represents the Llixha synclinal 
structure limited by the anticline structures of Elbasan-Valeshit 
in the east and Papri in the west. The Elbasan-Valeshi structure 
rises at the surface through the calcareous formation of Creta- 
Paleogenit, while in Papri the calcareous formations are at greater 
depth. The orientation of these structures is SW-NE as is the rest 
of the Albanides structure [2]. The formations in the region are 
mainly composed of flysch (three Oligocene and lower Miocene 
(less molasses) sections) while calcareous formations with lim- 
ited surface area are mainly composed of olystolit, through the 
over lined terygen rocks. The biggest part is composed of lower 
Oligocene (Pg3!) rocks, consisting of 6 lithological units as shown 
in Fig. 2 [1]. The lowest part of the formation is composed of thin- 
ner flysch units intercalated by clays and sandstone, while the 
upper part is composed of thicker flysch with conglomerates. The 
lithology of the middle and upper Oligocene formation is mainly 
flysch with differently shaped packs. Interesting is the lithologi- 
cal evidence of the upper packs of the middle Oligocene where 
among the flysches are located with lengths up to 100m. They 
are organized in a chain, preserving the general orientation of the 
other strata. Together with the conglomerate sandstone formation 
they form a chain alternating in relief. The lithology of the lower 
Miocene (Acquitan-Burdigalian) deposits is characterized by Mar- 
goles alternating with the clays-sandstones as well as the presence 
of sandstone-conglomerates. 

The Llixha syncline represents a depressed structure with 
eastern asymmetry filled in the central part with terygen, fly- 
sch and molasses deposits. The eastern part is distinguished 
by an easterly drop. The tectonics puts the upper deposits 
in contact with the lower flysch Oligocene formation and has 
made the surface intrusion of the Eastern anticline calcare- 
ous formations possible. This tectonic layout is regional and 
includes the Western anticline chain of the Kruja tectonic zone 
[3]. 

The Llixha system comprises a reservoir which feeds the 
southern part of the Shkumbini River. The main hydrological char- 
acteristic of the region is the presence of several hot springs. 
Their position is related to the Kruja geothermal zone and 
they are connected with the calcareous olystolits of the upper 
Oligocene conglomerate-sandstones. The temperature of the hot 
springs varies, and ranges from 50 to 68°C while flow rates vary 
from one spring to another, without any seasonal characteris- 
tics. 


2.2. The geological structure of Tirana 


Tirana aquifer is related with the syncline deposits of Tirana, 
whose morphology represents a depression of 10-12 km of width 
and 70-80 km of length, Fig. 4 (Geological map of Albania, 1984). 
The aquifer is composed of flysch, gravels, sandstones, alevrites and 
quaternary clays. Well yield in this region varies between 7 and 
101/s, while for the Tortonians molasses which are composed of 
sandstone, clays and alevrites the thickness is 60-150 m while the 
well yield is about 3-41/s. 
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Fig. 1. Map of Albania. 


3. Heat transfer theory of an elementary volume AV is equal to the thermal energy which 
crosses the surface S. A solid medium is considered, which is not 
3.1. The differential equation for heat transfer generating any energy (so the energy is only flowing through a 


surface S). The temperature T at a point P(x, y, z) will be a con- 
This equation is a mathematical expression of the first law of tinuous function of the position and time. For a homogenous solid 
thermodynamics, the energy preservation low. The heat increase medium, in which the thermal volume heat capacity is independent 
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Fig. 2. Geological map of the Llixha region. 


In polar coordinates: x=rcosa; y=rsina, hence: 


VT = 


T 139 187 T 1əT 


Ər2 ` rər r2 a2 ` 


az2.—s a. ott 


Let us assume that in point P(x, y, z) the energy input per unit of 
time and unit of volume is A(x, y, z, t), then Eq. (2) is rewritten as: 


A(x, y,Z,t) _ 1 0T 


27 4 
vi À ~ a ot 


(3) 


If thermal conductivity depends on position and temperature the 
equation can be written as: 


ar a (f.aT\ af.ar\ a(.ar 
PCE = Ox (az) y (5) E (a3) (4) 
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Fig. 3. Geological cross-section through the area (I-I in Fig. 2). 


Its solutions are often relatively simple and can be found for dif- 3.1.1. The basics hypothesis 
ferent functional forms of A(x, y, z). If the thermal properties are Let us describe mathematically the heat exchange process 
temperature dependent the equation is non-linear and numerical between flowing fluids in a wellbore with appropriate boundary 
methods can be used for its solution. formations. Let us assume: 


A-A profile Student City, Tirana 


Fig. 4. Geological map of Tirana aquifer. 


N. Kodhelaj / Renewable and Sustainable Energy Reviews 15 (2011) 2534-2544 2539 


e The formation is homogenous and isotropic; 

e The flow is in axial direction only; 

e The heat flow is only by conduction; 

e There is no radial temperature gradient in the wellbore; 

e In the boundary zone of the wellbore the formation temperature 
is a known function of depth; 

e The thermal properties of the fluid and formation are constant; 

e The radial temperature in wellbore is constant; 

e The yield from the wells is constant; 

e The flow in wellbore is vertical (1D); 

e The vertical heat conduction - is much less than the horizontal 
one, i.e. negligible. 


3.2. Initial and boundary conditions 


In order to solve the thermal diffusion equation let us assume 
that the formation temperature is function of the position (x, y, z) 
and time (t) [4]. In addition initial and boundary conditions must 
be specified. A specific moment of time is chosen as the origin of 
the time coordinate. At that time the temperature distribution is: 


T(x, y,Z,t =0)=f(x,y,z) or T(r,0,z,t=0)=f(r, 0,z) (5) 


If the radials symmetric case of a flowing wellbore is considered, 
as an example. After a certain production time the well is shut 
down. To determine the temperature distribution in the wellbore 
during the shutdown time, the end of production is considered as 
the origin of the time coordinate. In this case the temperature f(r, 
z) must be known, i.e. the initial conditions. To specify the tem- 
perature field of a medium the boundary conditions must also be 
known beforehand. Different kinds of boundary conditions can be 
specified: 


e Surface temperature is known. It can be constant or function of 
position and time, T=/{x, y, z, t). 
e The amount of energy flowing through the surface is known: 


qs(X, y, Z, t) = —À a { 2 = the derivation perpendicular to the surface } (6) 


Linear surface heat flow. In such a case the amount of energy 
transmitted through a given surface is proportional with the tem- 
perature difference between the surfaces and the surroundings. 


qs = &(Ts — To) (7) 


where T; is the surface temperature; Tp is the surrounding tem- 
perature; a is the heat transfer coefficient, and 


liMg—+ooTs = To (8) 


e The connecting surface between two media with conductivities 
A, and Az, respectively. If T4 and Tz are the temperatures of the 
media, then 


Tils = Tals; M 


OT. OT. 
7 Is (9) 


2 
an 8 = 42 gn 
3.3. Temperature field modeling 


In order to create a physical model to enable the stable state of 
heat flow through a formation to be studied, the analogy princi- 
ple may be used. Let us assume that n dimensionless parameters 
Ij,...J, describe the heat conditions in the formation and sur- 
roundings. In such a case, the model results should be expressed 
as n dimensionless parameters I} , ..., I}. The conditions that must 
be fulfilled in order to use the model in a real reservoir are I; = 
D,- h =I). For stable heat transmission, the analogy with the 
electrostatic field can be used. For the 2D case and with the thermal 


conductivity À = A(x, y), the boundary conditions are: 


T(z1,x)= Ti; T(zz,x)=T2; -l<x<l; 


_ OT(Zz,x), 
= Se =0 (10) 


In Table 1, a summary of the analogies between hydrodynamic, 
heat transfer, electrostatic and electricity transmissions fields is 
presented. As we can see from Table 1, the temperature change 
AT is analogous with the potential difference AU. Assuming the 
length unit, /;, the dimensionless coordinates are: 


l= (11) 


The relationship à =A(x, z) is analogous with p= (x, z). A partic- 
ular care should be used in unstable field modeling. Considering 
the hydrodynamics analogy with the thermal one, in particular 
the transitory flows of an incompressible fluid through a porous 
medium, the main characteristics of this regime are the high values 
of the hydraulic diffusivity [2]: 


k 


= ona (12) 


n 


Its thermal field analogue is the thermal diffusion coefficient (a). 
The analogies principle can only be used in the case when the 
dimensionless times are equal. In thermal field modeling, the fact 
that the temperature is not affected by surface topography and 
groundwater flow is very important. To determine the boundary 
conditions temperature maps at depth are widely used [5]. 


3.4. The unstable temperature field 

Before starting investment in a geothermal project, the stability 
of the temperature field involved, in space and time, needs to be 
confirmed as well as the projects overall sustainability. The meth- 
ods that can be used to answer such questions are numerous, but 


here we will apply the finite element method [6]. Considering the 
differential equation: 


ð (p oT _ 9 p, 2T „2 (p T ee 
Ge oe) ay kOe oe ee | ee 
The functional of this equation is [6]: 
2 2 2 
_ _ ky [ oT ky / oT kz [ oT 
ramon S S SE) -3 s(a) 
oT 
xarcyae+ | | (a-:3) dxdydz (14) 
m ot 


Let us divide this functional into two parts: 
2 2 2 
2 ky / OT ky ( oT _ kz oT . 
w= SSE z(a) +2 (Gz) |e 
Xr = iy (0-<F) dxdydz (15) 
2 ot 


To integrate the Eq. (14) let us assume the time interval (t, t+ At) 
through into dT/dt=C 


0 
(13) 


co = 


2540 


Table 1 
The analogy between different physical fields (based on [4]). 
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Hydrodynamics Heat transfer 


Electrostatics Electricity transmission 


Pressure Temperature 

P T 

Pressure gradient Temperature gradient 
—yP -vVT 


Permeability/viscosity k/u Thermal conductivity 4 


Heat amount transfer q =—VT 
Isothermal surface T= C 
Split surface {Æ = 0 


Velocity vector v=-k Vp 
Isobaric surface p=C 
Impermeable layer w =0 


Electrostatic potential Potential 

@ U 

Field vector Potential gradient 
E=-v@ -vU 

Dielectric constant ¢/47 Specific resistance p 
Dielectric displ. (e/4n)E = —(€/47)V® Current | = -VU 


Surface with electrostatic potential ® =C 
Force line 22 = 0 


Surface with potential s =C 
Potential line 3% = 0 


T;: temperature values in the nodes [=1, 4 


Ty 
OX OXy Xr. 1 T 
IT, ont OT T= 5A { Ny Nz N3 N4 } T where 
T4 
N(x, y, z) = (a; + bix + ciy + diz); i=1,4 (16) 


The partial derivates of the functional can be calculated as follows: 


ox 

a b? bı bz bı b3 bı b4 A 
ox E OT> _ kx | boby b2 bəb} bb ky C2C1 
ƏT; F oX 364 b3 by b3 b2 b2 b3 b4 364 C3C1 

aTa b4b1 b4b2 b4b3 b3 C4C1 

ox 

OTs 


T4 

I / f N?dxdydz 
2 

J | | senavave 
2 

J | | sonravave 
2 

J | | sesavave 
2 


Let us make the following replacements: 


b? bib} bib} by bg ci C&C C3 C4 

kx bob; b2 bəb3 bzb4 ky C2C4 c2 C2C3 C2C4 

364 b3bı b3b2 b3 b3b4 36A C3Cy, GG ce C3C4 
bab, b4b2 b4b3 b4 C4C1 C4C2 C403 c3 


QA 
<= = [F] 
where [H] is the matrix of conductivity of the thermal field, [P] is 
the matrix of instability of the thermal field, [F] is the source vector. 
After these replacements the definitive form of Eq. (17) 


becomes: 


(20) 


[H] |T| Heri- |F(t)| =0 (21) 


J | | NNadxayez J | | NNsaxayaz 

2 Q 

J| | pacavae J | | MNsdxayaz 
2 Q 

J | | senatvavac J| | Miye 
2 Q 

J | | sesetvavae J | | NNsdxayaz 
2 Q 


J | | Miye J | | simatvavae J | | simsavavee J | | NNadxiyaz 

Q Q R2 R 

J | | NNidxiyaz J| | Mixa J | | MNsdxiyaz J | | sematvava 
Q Q 2 Q 

J | | sonavave J | | sometvavae J| | Miye J | | somatvava 
Q Q 2 Q 

J | | NNiaxayaz J | | NNadxayez J | | NNsaxayaz J| | Niixayaz 
R2 R R2 R 


4. District heating system using the geothermal water 


In general the geothermal water for district heating systems is 
taken directly from low temperature reservoirs. Another way is to 
use the geothermal water through the heat exchangers to heat up 
the fresh water. The hot water can be stored in tanks if appropriate. 
This water then is transmitted to the buildings and can be used for 
heating and as tap water. Heat flow to the buildings is controlled by 


Cic C13 CG d? dıdz dıd3 dıd4 Tı 
c2 GCs GG kz dod ds dzd dod, Tə 
C3C2 c C3C4 36A d3dı d3d> d d3d4 T3 
C4C2 C403 A dadı d4dz d4d3 d T4 
/ | | N,Nadxdydz aT, 
A “Ot 
I i f N>N4dxdydz | | 3T2 1 0 
Q ar | QAJ1| jo 
IT 4 /|1| Jo 
J | | sonatvava Or 1 (0) 
R 
OT4 
J| | Mitya “Ot 
2 
(17) 
d  dydy did; dıd4 
k, dod, d2 d2d3 dz2d4 
364 |d3d; d3d2 a d3d4 [H] (18) 
dad; da4dz d4d3 d3 
=[P] (19) 


the mass flow. In the following are given the basics calculations for 
the geothermal district heating. The main elements of the geother- 
mal district heating are the radiators and also this use is affected 
by the water thermal energy, building heat loss, pipe heat loss and 
building energy storage. A short description for all of them is given 
below. 
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Table 2 
Estimated water temperature at depth based on different geo-thermometers. 


Geothermometer Spring “Nosi” Llixha, Elbasan (°C) 
Fournier 254 
Truesdell 235 
Na+Ka+Ca 143 
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CROSS SECTION 
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Fig. 5. Geomorphologic map of the Llixha region. 
4.1. Calculation procedure for house heating systems 


Radiators are the heat exchangers and they make the trans- 
fer of the heat from the geothermal water to the indoor possible. 
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Fig. 6. The finite volume grid. 
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Fig. 7. The temperature magnitude. 
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Fig. 8. The density magnitude. 


The relative heat capacity of the radiator is given by the following 
relation: 


Qrad AT in 4/3 

ge ~ (a 

where the subscript 0 denotes the design conditions. The temper- 

ature difference ATm is calculated [7]: 
(Ts-T)-(r-T) _ Ts-T 

In ((Ts — T;)/(Tr -T;)) 1n ((Ts — T;)/(T; — T;)) 


The water thermal energy which the geothermal water will give to 
the indoors, through the radiators is: 


Qrad = cpm(Ts — Tr) (24) 


(22) 


ATm = (23) 


The relative thermal energy given by the geothermal water is [7]: 


om _ _m(Ts ~ Tr) 
Qa  mo(Tso — Tro) 


(25) 
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Fig. 9. The velocity magnitude. 
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Table 3 
Thermal characteristics of the rocks. 


Age Lithology Thermal conductivity (W/m °C) Thermal resistance (m°C/W) 
Tortonian Sandstone 1.6 0.625 
— Clay LS 0.666 

Table 4 


Physical characteristic of the sandstones. 


Layer Filtration coefficient Water conductivity (m/day) Specific yields (Im/s) 
Tortonian sandstone, upper part 0.08 3.2 0.038 
Tortonian sandstone, lower part 0.08 5:3 0.038 


The biggest part of the energy loss in such system may come from 
the buildings. These losses can be calculated as [7]: 


Quoss = ki(T; a To) (26) 


where building heat loss factor k; is constant. The relative losses are 
[7]: 

Qioss — Ti = To 
Qiosso Tio — Too 
The other losses in such system are the pipe heat losses. To deter- 
mine the amount of losses in the pipes it is necessary to know the 


pipe transmission effectiveness parameter T. Its value [17] is given 
by: 


Ts — Tg Up 
T= T exp ( me (28) 


The reference value of t can be calculated from the reference flow 
conditions as follows: 


Tso — Tgo Up 
To = = ex 29 
0 rs p moc (29) 
The values of Up and cp are assumed to be constant in the system. 


By combining Eqs. (28) and (29), the transmission effectiveness can 
be calculated: 


(27) 


T= grom (30) 


The supply temperature to the house can be calculated: 


Ts = Tg + (Tı — Tg)t = Tg + (Tı — Tg )t{o/™ (31) 
and the return water temperature at the pumping station is calcu- 


lated: 


T2 = Tg + (T; — Tg)t = Tg + (T; — Tero (32) 
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Fig. 10. The velocity vectors. 


A building’s heat storage (dependent on the amount of the thermo- 
insulation) can be very helpful for the heating system. Normally the 
builds are like a heat store. The building heat storage is calculated 
by [7]: 
dT; 1 
dt c 


Qnet = = (Qsup — oss) = [MepCTs — Tr) — ki(Ti — To) 
(33) 


5. Borehole heat exchangers 
5.1. Equation for heat flow 


The equation for rate of heat transfer from the fluid in the heat 
exchanger to the earth mass is: 


- =UxAT (34) 
where Q is the rate of heat exchanger (BTU/h or W) for the whole 
heat exchanger length; L is the length of heat exchanger (m); U is 
the conductance rate for heat transfer from the circulating fluid to 
the earth (BTU/h/°F or W/°C/m) for the operating conditions: AT is 
the difference in fluid temperature: 

= T+T 


AT=2 1- (35) 


Tp is the earth temperature (°F or°C); T4 is the fluid entry temper- 
ature (°F or°C); Tz is the fluid exit temperature (°F or °C). 
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Fig. 11. Calculations for the radiators. 
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Fig. 12. Lithological column of a well drilled in Tirana. 


The conductance term U for heat flow from fluid in the heat 
exchanger to the earth can be estimated with the conductance 
coefficient for composite cylinders. The impedance to heat flow 
is caused by the thermal resistance of the pipe wall and the soil 
cylinder around the casing. Fluid surface resistance films are small 
relative to the other terms and are encompassed in the two resis- 
tance terms. We can express U as: 


27 27 
-~ soil resistance + pipe resistance Rs + Rp 


(36) 


6. Results of simple hot spring modeling 


The six hot water springs at Llixha in the Elbasan region, have 
water temperature up to 65°C and flow rate up to 231/s. Esti- 
mated reservoir temperature at depth associated with the hot 
springs, based on the chemical composition using different geo- 
thermometers, is given in Table 2. 

These values show that the water is coming from great depth 
where the average temperature is over 200°C. As it is shown in 
Fig. 5 [8], the hot springs are situated in the middle of the village. The 
hot water has been used only for balneology for several centuries 
possibly since the time of the Roman Empire. The first modern use 
dates back to 1937 with the building of the “Hotel Park” medical 
centre. The use of the water flowing from these springs can help 
to improve the economical effectiveness district heating in the vil- 
lage. The realization of such a project will allow the utilization of 
geothermal water as an energy resource for the first time in Alba- 
nia. The purpose of the calculations presented in this chapter is to 


show that the water of these hot springs can be used for the district 
heating of the village community. 


6.1. The 3-D modeling of the temperature field of Llixha, Elbasan 
region 


A finite volume model was set up for a crustal volume with an 
area of 10km x 10km and 5km thickness to model the temper- 
ature, density and fluid velocity distribution in the Llixha region. 
The grid is shown in Fig. 6. Here it is assumed that the medium 
is homogeneous and isotropic and that ky=ky=k,=2W/mK [9]. 
We also know that Q=201/s (corresponding to m;=Q/6=3.31/s or 
3.224kg/s for each of the hot springs), cp =4180J/kg°C [8]. The 
temperature at depth in the formation is set at 221°C while the 
temperature of the water at the surface is in the range 60-65 °C. 
The temperature gradient of the surroundings is assumed 12 °C/km. 
The modeling software FLUENT is used to solve the problem, it pro- 
vides calculation results for temperature, density and velocity for 
the volume modeled. In the model water flows with a velocity of 
1.25 x 10-7 m/s. The results for temperature, density and velocity, 
as well as velocity vectors, are shown in Figs. 7-10. 


7. Basic district heating design 
7.1. Calculation for the radiators 


The steps used to calculate the radiator parameters are given 
in Section 4. In these calculations, it was assumed that the sup- 
ply temperature of the water is 60°C, the reference temperature 
65°C, and the ground temperature 6 °C. The calculations were done 
for 4 different scenarios: The indoor temperature is assumed to 
change at the range of 18-20°C, the outdoor temperature in the 
range of —10 to —4°C, the return water temperature in the range of 
33-40°C, the reference inflow in the range of 3.93-7 kg/s, the ref- 
erence system inflow 18-24 kg/s. Based on these data, the relative 
heating of the radiators was calculated as 0.83-0.92, the relative 
heating of the building as 0.87-0.88, and the transmissivity coeffi- 
cient t = 0.94. To be within these parameters, it is sufficient that the 
supply water temperature be 60°C, and the inflow of the system 
16-27 kg/s. Thus the parameters for the Llixha thermal springs sat- 
isfy all these demands. Fig. 11 shows the results of the calculations 
for the radiators supply and return water temperatures vs. indoor 
temperature. 


7.2. Borehole heat exchanger calculations 


Heat exchanger length is calculated through Eq. (34). The ther- 
mal characteristics of the rocks drilled by the well of Fig. 12 are 
given in Table 3, while the physical characteristics of Tirana torto- 
nian sandstones are given in Table 4. 

For this well is known: 


e Dynamic level 100 m; 

e Static level 50m; 

e Level falls 40m; 

e Pump depth 108 m; 

e Average yield 11/s; 

e pH 7.5-8; 

e Water temperature 13-15°C; 

e General strength 10-15 °Gj; 

e Mineralization 250-350 mg/l; 
e Formula HCO3-Ca-Mg; 


The thermal resistance for the formations Rp drilled by the well 
is calculated as follows: 


(40 x 0.666 + 34 x 0.625) _~0.647™ xC 


Rp = 40+35 W 
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Fig. 13. The economical analyses of three different systems. 


The conductance rate U is calculated: 


20 27 W 
u= Rp 0.647 =a) Nl eae 


According to the measurements, the earth temperature in 100m 
depth in Tirana is Tg =18°C. The fluid exit temperature is 3.5°C 
lower than the fluid entry temperature. The difference temperature 
of the system is calculated: 


h+h _6§75°¢ 


AT = To — 
For these parameters the heat exchangers, for the installed capacity 
of 100 kW result: 


Q 100 x 10° 
“Ux AT 9.711x6.75 


So the problem is solved through drilling process of 15 wells 100 m 
depth each. The economic analyses (profit-expenditure) of three 
systems: Water-Water geothermal pumps, Earth-Water geother- 
mal pump and the existing coal heater are presented in Fig. 13. So it 
can be clearly seen that the maximal profit is for the Water-Water 
system than for the Earth-Water and finally for the coal heater. The 
payback period for the geothermal heating pumps varies from 2.8 
up to 10.75 years. 


L 


= 1525m 


8. Conclusions 


Based on the calculations presented in this report the following 
can be concluded regarding utilization of the hot springs of the 
Llixha-Elbasan hot-spring area in Albania: 


e The water temperature is expected to be stable in the future; 

e The geothermal water from the Llixha hot springs fulfils all 
requirements for district heating’s in the region; 

e Considerably higher temperature expected through further well 
drilling; 

e Albanian geothermal regime allows different scale borehole heat 
exchangers applications; 

e Demographic and geological features of the student’s city allow, 
and furthermore are feasible the borehole heat exchanger’s uti- 
lization. 


References 


[1] Kogiaj S. On the construction of the Earth Crust in Albania according to the first 

on set of “P” waves in the seismologic station. Bull Geol Sci 1989;1:400-10. 

[2] Aliaj S, Hyseni A. The neotectonic map of Albania, scale 1:200000. Tirana: Sh. 

B. L. U.; 1996, 87 pp. 

[3] Hyseni A, Melo V. The geodynamics of new movements in Albania and their 

influence in resources and environment. National programmer for research 

development, Toena, Tirana; 2000. 

[4] Plummer C. Physical geology. N.Y.: Elsevier; 1988, 580 pp. 

[5] Cermak V, Haenel R. International association for seismology and physics of the 

earth’s interior; 1988. p. 30-3. 

[6] Osmani S. Numerical methods. Tirana: Sh. B. L. U.; 1997, 125 pp. 

[7] Nappa M. District heating modeling. University of Iceland, report; 2000. 47 pp. 

[8] Frashëri A, Cermak V. The geothermal atlas of Albania. Tirana: Sh. B. L. U.; 2004. 

p. 65-89. 

[9] Osmani S. Finite elements, applications I. Tirana: Sh. B. L. U.; 1987, 210 pp. 

[17] Valdimarsson P. Modelling of Geothermal District Heating Systems, Ph.D. The- 

sis, University of Iceland, Faculty of Engineering, Reykjavik, Iceland 1993, 121 


p. 


Further reading 


[10] Buffon GL. Histoire naturelle, générale et particulière - 651 S. De l'imprimerie 

royale, Paris. Commented edition: Roger J. (1962): Buffon, Les époques de la 

Nature - 495 S. Mémoires du Muséum National d'Histoire Naturelle, Série C, 

Tome X, Editions du Muséum, Paris; 1778. 

[11] Harlow JH, Klapper GE. Residential heat pump experiments in Philadelphia- 

installation and operating experience. AIEE Trans 71/II, New York; 1952. p. 

366-75. 

[12] Ingersoll LR, Adler FT, Plass HJ, Ingersoll AC. Theory of earth heat exchangers 

for the heat pump. ASHVE Trans 56; 1950. p. 167-88. 

[13] Osmani S, Korini Th. Geostatistics. Tirana: Sh. B. L. U.; 2001, 310 pp. 

[14] Turcotte D, Schubert G. Geodynamics. Cambridge; 2003. p. 279-87 (580 pp.). 

[15] VDI. Thermal use of the underground - ground source heat pumps. Guideline 

VDI 4640, part 2. Berlin: Beuth Verlag; 2001. 

[16] WuY, Pruess K. An analytical solution for wellbore heat transmission in layered 
formations. SPE 17947; 1988. p. 350-500. 


